cpi uiaii 



Toward a General Understanding of the Scaling Laws in Human 
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Abstract. - Recent research highlighted the scaling property of human and animal mobility. An 
interesting issue is that the exponents of scaling law for animals and human in different situations 
are quite different. This paper proposes a general optimization model, a random walker following 
scaling laws (whose traveling distances in each step obey a power law distribution with exponent 
a) tries to diversify its visiting places under a given total traveling distance with a home return 
probability. The results show that different optimal exponents in between f and 2 can emerge 
naturally. Therefore, the scaling property of human and animal mobility can be understood in 
our framework where the discrepancy of the scaling law exponents is due to the home return 
constraint under the maximization of the visiting places diversity. 
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The scaling laws of animal mobility indicate a class of 
random walks (which is also called Levy flights), whose 
each step jump distance I typically follows a power-law 
distribution of 



p{i) oc r 



(1) 



where 1 < a < 3. This kind of random walks capture the 
property of the foraging behaviors [?] of albatrosses [I], 
terrestrial animals [2j and submarine predators [3]. On 
one hand, the exponents in animal foraging behaviors are 
always approaching 2 (e.g. a « 2 for albatross PQ, deer 
[2], and Atlantic cod (Gadus morhua) [3], and, a ~ 1.9 
for lcatherback turtle (Dermochelys coriacea), a ~ 1.7 for 
Magellanic penguin (Spheniscus magellanicus)). Even for 
human hunters, as Brown et al. [A pointed out, the scaling 
exponent is also approaching 2 as other foraging species. 

On the other hand, as recent experimental researches 
highlighted, the scaling laws (EqUJ can be generalized as 
an approximation to the human traveling patterns |TJ|5H7] . 
However, the exponents are much less than the ones of an- 
imal foraging behaviors, such as a s=s 1.59 in Brockmann's 
bank notes tracking study [5], a ~ 1.2 in Gonzalez's mo- 
bile phone users mobility [6] and a ~ 1.55 in Song's high 
resolution data [7J. 

It is notable that, the difference of scaling exponents 
is not negligible (Acti/an « 25.78% for individual or 
Aa p /a p w 66.6% for population [6]) and is mainly due to 



the behaviors in different situations. For the human trav- 
eling or transportation behaviors (e.g. bank notes and mo- 
bile phone users), the exponents are approaching 1; while 
for the active searching (foraging) in the food lacking en- 
vironments 019], the exponents are around 2. 

What is the mechanism underlying the scaling behaviors 
with different exponents? On one hand, Visawanathan et 
al. have proposed a model [10] to explain the animal's 
foraging behavior, in which they assumed that the forager 
searches for food "sites" following EqQ] with variant ex- 
ponents. They calculated the efficiency of search as the 
mean flights taken between two successively visited sites 
and obtained the optimal exponent a op t — 2 which yields 
the best searching efficiency. Although this model can fit 
the foraging behaviors very well, it fails to explain human 
traveling behaviors (a ~ 1). 

On the other hand, some recent studies tried to explain 
the human traveling patterns, i.e., the scaling behaviors 
with exponent approaching 1. For example, Song et al. [7J 
assumed that in each time step the individual may look for 
a new place which is never visited before or go back some 
visited places (e.g. home or work place). This stochastic 
model can explain their observed human mobility data. 
Han et al. [11] attributed the origin of human mobility pat- 
terns to the hierarchical structure of streets. Nevertheless, 
these models cannot explain the foraging behaviors. 

Therefore, it still requires a general explanation. De- 
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spite the significance of the discrepancy of scaling expo- 
nents is pointed out by [T2], it is still poorly understood. 
In this paper, we try to propose a general model of human 
and animal mobility. 

First, we assume that one of the most important driv- 
ing forces of universal mobility patterns is diversity. That 
is, the random walker "tries to" maximize the diversity of 
his visiting locations (which is measured by the Shannon 
entropy) within the limited total distance [13] . This con- 
straint can be also understood as the energy constraints 
of human traveling behaviour [?. ?.I14]. Here, we should 
claim that the maximum diversity (entropy) is achieved 
not by the intentional calculations of the random walker 
but as a most possible result of the underlying stochas- 
tic process [13, 15, 16 . Furthermore, this assumption is 
consistent with Visawanathan et al. 's model in the forag- 
ing behaviors because seeking foods with highest efficiency 
in the food lacking environment is equivalent to diversi- 
fying visiting places under the total distance constraints 
[10] . And also, diversity is a very important driving force 
in people's daily lives and economic development [17H19) . 
This can be reflected by maximizing the diversity of vis- 
iting places in human traveling. This assumption is also 
coincident with exploring new places in Song's model [TJ. 

Second, the other important driving force exerting on 
the random walker is the home return constraint. In other 
words, the random walker must return home (a giving site) 
with a fixed probability in our model. This assumption is 
also supported by our daily experience and previous works 
[TJ. Both animals and human beings always return some 
fixed points (homes, working places, schools etc.). 

Under the above two important presumptions, we ob- 
serve that different scaling exponents in between 1 and 2 
can emerge naturally under different home return proba- 
bility in the process of maximizing the diversity of visiting 
places. Therefore, our model provides a possible general 
understanding of human and animal mobility scaling prop- 
erties. 

The Model. — Suppose the mobility space is an L x L 
toroidal lattice (where L is a given constant denoting the 
world size). The random walker indicating an individual 
of human or animal can travel around in this lattice ac- 
cording to the given scaling behavior. In other words, the 
movement pi = {x, y) in the ith time step is a two dimen- 
sional random vector whose length \p t \ following Eq. [1] 
with the exponent a. The purpose of this model is not 
to explain the origin of scaling law but to study how the 
home-return tendency affects the power law exponent a, 
therefore, scaling law, i.e., Eq. [T] is a basic assumption 
which can be supported by the previous studies [51[6lll0j. 

Second, we should consider the home return constraint 
which is a distinct feature of our model. In reality, human 
or animal individual periodically visits or returns some 
specific places such as home or working place [7J. How- 
ever, in our model, we suppose that the random walker 
should return to only one place called "home site" with a 



fixed return probability r. The reasons why we select only 
one site as home are following: 1. To keep the model as 
simple as possible (The multiple fixed sites cases will be 
studied in the future); 2. In the case of human traveling, 
although people always travel in between home and work- 
ing place, the distance between the two points is quite 
small considering the whole mobility space 5 , 7 . So, we 
can use "home site" to stand for the aggregation of these 
two points. Without losing generality, we suppose the 
"home site" is just the axis origin point O = (0, 0). 

Thus, the integrate stochastic process can be described 
as follows. In time step t, suppose the current position of 
the random walker is x t . Then the next position in the 
time step t + 1 is 

J Xt+pi, if xt = O or 9>r (2) 
Xt+1 ~ \ O, if t = or act ¥= O and 6 ^ r.(2') 

where, 9 is an independent random number evenly dis- 
tributed in [0, 1] and pi is an independent random vector 
in two dimensions with a random length distributed fol- 
lowing the scaling law (Eq. [TJ. So the random walker 
starts from the origin at t = and jumps out to other 
places with probability 1 — r following Eq. [T]step by step. 
It may return back home with a probability r in each time 
step. 

To avoid the infinite repetition of this model, a stopping 
condition is added. We assume that the total distance of 
this random walker instead of the total number of time 
steps is given. That is, we have the following constraint 

T 

J2\pt\<W = cL, (3) 
t=i 

where W is proportional to L (c is a given constant). It 
is reasonable to assume that the total distance is propor- 
tional to the total energy consumed by the random walker, 
so this constraint is consistent with the energy consump- 
tions in spatial network [?, ?,[14l[20]. Actually, for given 
W, the total time steps T is a random number. Theoret- 
ically, the distribution and expectation value of T can be 
obtained by solving the first passage time problem, and 
in this way the average stopping time T is dependent on 
a(see the supplementary material). However, in the sim- 
ulation, we just generate a set of random vectors p t until 
the total distance constraint Eq. [3] is violated to get T. 

Eq. [TJ [5] and |3] actually define the dynamical pro- 
cess of the random walker. Under the dynamics, dif- 
ferent sites may have different chances to be visited by 
the random walker in his whole life. Obviously, the sites 
around the home site will be visited more frequently in 
the tight home return constraint (i.e., large r). Formally, 
let {x~i, x~2, x~t} denote the sequence of sites that the 
random walker have visited in all steps of one realization 
of the random process. Then, the visited probability of 
each site x is defined as: 

T 

pW=Y,S*-3/T, (4) 
t=i 
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where, 6g is the Kronccker's delta function. It is 1 only 
when y = 0, otherwise it is always 0. We can plot visited 
probability of all sites in four simulations in Fig [TJ This 
figure can be compared to the empirical distribution in [B] . 
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Fig. 1: The visited probability of each site in four different 
cases with L = f 00, W = 1 00000. The visited probability dis- 
tribution is comparable with the empirical results in [6j quali- 
tatively 

Finally, we assume that the underlying reason to form 
mobility scaling patterns is the driving force of maximiza- 
tion of the visiting places diversity. This trend can be 
quantified by maximizing the Shannon entropy of the vis- 
ited probability. The Shannon entropy is defined as |21) : 



S 



LxL 

£ 

X 



p{x)\ogp(x) 



(5) 



The summation is taken for all sites. 

Therefore, a complete optimization problem is defined. 
We will find an optimal exponent a to maximize Eq. [S] un- 
der the dynamical rule Eq. [2] and total distance constraint 
Eq. [3j Obviously, different home-return probability r will 
systematically influence the optimal exponent a. Next, 
we will show the simulation results of this model in the 
main text. And the analytic results in some special cases 
are presented in the supplementary material. 

Results. — Because our main purpose is to study how 
the home-return pattern affects the optimal exponent. We 
can first consider two extremal cases, i.e., "never back" case 
(r = 0) and "immediate back" case (r = I). In the first 
case, the trajectories of the random walker are not con- 
straint by return probability. It means that the home site 
is not effective and the random walker will travel until 
W was fully consumed. In the second case, the random 
walker immediately flies back home after one step walk. 
The optimal exponents maximizing the entropy should be 
quite different in the two cases. The simulation results are 
shown in Fig [5] 

From FigfJl we can observe that the optimal exponent 
Oiib in the "immediate back" case (r — l,W = L,10L) 
is around 1. However, the optimal exponents a n b for the 
"never back" case (r = 0, W = L, 10L) are from 2 to 2.3 
with different W . The optimal exponent in other cases of 




Fig. 2: Entropy S v.s. exponent a in extremal cases (L = 
10000, IF = L or 10L). Each result on the figure is the en- 
semble average of 100 experiments. We can observe that the 
information entropy of visited probability will change with a. 
The optimal a which can maximize S appears in different po- 
sitions with different home-return probability r. The optimal 
exponent is around 1 when r = 1 and 2 when r — 0. 



the optimized problem should be in between of ctib and 
&nb- With simulation results in the above two extremal 
cases alone, we have confirmed that the home return prob- 
ability r does have influenced the scaling law. 

To approach how the home-return probability affects 
the scaling law, we have done another group of simulations 
in which the overall distance W = 10 * 10000 is fixed and 
home return probability r varies ranging from to 1. We 
draw the curves showing how the exponent a affects the 
entropy S for different r in FigJ3]and|31 

According to Fig. [3]and|H the optimal exponent a opt 
changes from around 1 to 2 as r decreases. This result tells 
us that to maximize the total entropy, a random walker 
seldom travel long distances when it is more bound to the 
home site. And as we expected, the information entropy S 
decreases when r increases, since random walker's capacity 
of information foraging is constraint by home-return. 

One may ask a question: how the home-return proba- 
bility (i.e., r) affects the optimal exponent a opt ? To study 
this issue, we conducted another group of simulations with 
different r and plot the curve of (r, a opt ) in Fig. [5j From 
the plot we see that the value of the optimal exponent 
a opt decreases slowly at first when r is in [10~ 5 , 10~ 4 ] and 
quickly falls down from 2 to 1.5 when r is in [10 -3 4 , 10 -1 ]. 
At last, it drops from 1.5 to 1 slowly when r is larger 
than 0.1. This result indicates that the optimal scaling 
exponent a op t is determined by the home-return proba- 
bility, especially when r £ [10~ 3 ' 4 , 10 -1 ]). It is reason- 
able to regard the home-return probability r — 0.1 as the 
critical point to distinguish foraging behaviors and travel- 
ing behaviors. The interval (0,0.1] corresponds to active 
searching behaviors and (0.1, 1] is for human traveling be- 
haviors. This critical point can be also observed in the 
plot of S op t (the corresponding entropy in the optimal ex- 
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Fig. 3: The dependence of entropy S on a in different r(L — 
10000, W = WL). The entropy S changes with a and get 
maximal values in the interval around [1,2] when r changes 
from to 1. All curves are average results of 100 experiments 




Fig. 4: The dependence of entropy S on r and a (L = 
10000, W — 10L). All curves are average results of 100 ex- 
periments 
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Fig. 5: The dependence of optimal a and S on r. Where, 
L = 10000, W — L, 10L. We can observe that there is an 
apparent transition of the optimal a from 2 to 1.5 when r 
is around [10 _ , 0.1]. This transition is just the key factor 
that distinguishes human traveling and foraging behaviors. All 
curves are average results of 100 experiments 



ponent a) versus r. Notice that the S op t-r curve has a 
sudden change when r < 0.1. Furthermore, from Fig. [5] 
we can also observe that the optimal exponent can be dis- 
tinguished much easier in different r values when W is 
big. So, our conclusion is much better for the asymptotic 
behavior when W goes infinity. 

Conclusion. — To summarize, we propose a general 
framework to understand the mobility of human and an- 
imals based on the maximization of information entropy 
and home-return constraint. Our model reveals that the 
observed scaling exponents (1 for human traveling and 2 
for foraging) are the results of maximizing information en- 
tropy. Second, our model points out that the home-return 
constraint is very important and influences the optimal 
exponents of mobility scaling law dramatically. The dif- 
ference of scaling exponents between human traveling be- 
haviors and foraging behaviors is due to the strength of 
home- return constraint. 

As we know, human or animal mobility pattern is a 
key factor to study various related topics, such as spatial 
networks . 2"2~] gene pool diversity, mobile viruses and epi- 
demic spreading. Although we have not seriously fitted 
the simulation to empirical data, our exploratory work 
contributes to mobility research by highlighting the im- 
portance of home-return pattern based on a reasonable 
defined model and large scale simulation. We inferred that 
the home-return pattern is also crucial for approaching the 
origin of patterns for cities, traffic systems and other re- 
lated area. 

* * * 

We wish to thank Prof. Shlomo Havlin for some useful 
discussions. The work is partially supported by NSFC 
under Grant No. 70771011,61004107 and 60974085. 

Appendix. — In the main text, we have proposed an 
optimal model to explain the discrepancy of the power law 
exponent in human or animal mobility patterns. However, 
we only presented the simulation results due to the com- 
plexity of the model. To understand the mathematical 
essence of the model and get more solid conclusions, we 
try to analyze this problem in a more mathematical way. 

You will see, in some extreme cases, the model can be 
described by pure mathematical equations. Unfortunately, 
not all of these equations are solvable, so we have to give 
the numeric results instead of the exact analytic solutions. 
The main purpose of this analysis is trying to give a much 
clearer understanding toward the original problem but not 
a complete mathematical solution. We will discuss this 
problem in two extreme cases according to the home re- 
turn probability r, namely, r = 1, r = 0. 

Immediate Return (r = I). According to our dynam- 
ical rules, when r = 1 the random walker will jump out 
following power law distribution (Eq. 1), and return back 
to the home site immediately, and then, it will jump out 
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again and fly back, and so on until the T th time step when 
its total distance is consumed up. 

Because each jump-out and return back cycle is inde- 
pendent on the previous cycle, the final result can be con- 
verted to an equivalent problem: T random walkers start 
to jump one step out from the home site simultaneously 
according to the scaling law (Eq. 1). Therefore, the orig- 
inal temporal experiment is converted to an ensemble ex- 
periment. 

Suppose the probability of the random walker visits any 
site x = (x,y) is p(x), and the site x is visited by 77(2?) 
random walkers in the all T walkers. Therefore, the prob- 
ability that r)(x) — i is a binomial distribution, 



P{r,(x)=i} = C* T (p(x)y(l-p(x)) 
Thus, the average visiting times of site x is: 

(77(f)) = Tp(x), 
So, the average visit frequency of site x is: 

(r)(x)) 



T-i 



H{x) 



p{$), 



(6) 



(7) 



(8) 



which is independent on the total number of walkers T. 
Then, the Shannon entropy can be calculated as: 



LxL 



S = - ^ P(&) l °SP(x) 



(9) 



and, we have, 

dS(a) _log^ a+l + Ca-lJlog^- 



da 



a — 1 



(a -If 



(13) 



which is always smaller than when a > 1. Therefore, 
we can conclude that S(a) is a monotonic decreasing 
function. 

ii. When < a < 1 

We know that the distribution p(x) is not a standard 
Pareto power law distribution but have an upper bound of 
\x\ in this case, so the entropy cannot be calculated as the 
previous case. We will analyze the asymptotic behavior 
when L — > 00. 

At first, we know that p{x) is proportional to | a?| ^ 1 , 
therefore, Eq. [10] is still hold. However, \x\ cannot go to 
the infinity but have an upper bound L. And also, Z is a 
normalization constant, it is calculated as: 



L>\S\>1 

2ir 



2tt pL 

d6 
Ji 



\x\-° t - 1 dx'= I d9 I p~ a ~ x pdp 



(14) 



a-1 



(1-L 1 



So, when L — > 00, because a < 1, we know that: 



which is also independent on T. Next, we will give 
the concrete mathematical form of p(x) so that the 
relationship between S and the exponent a will be given. 
We know that p{x) is a Pareto power law distribution 
only when a > 1, L — > 00, so, we have: 
i.When a > 1, 



p(x) = \ \A 

for any \x\ ~^ 1, where 



1 



1 > 2 + y 2 )-^, 



(10) 



Z = I I \.r 

I\S\>1 

2tt 



•2tt z'+oo 

-<- 1 dx= I dO I p~ a ~ x pdp 

h (n) 



a-1 



Notice that there are 2tt\x\ points having the distance |x| 
from the origin and having the visit probability propor- 
tional to \x\~ a , so for each point x, the visit probability 
should be proportional to |a?| _Q_1 . 

Thus the entropy can be approximated by an integration 
when L — > 00: 

S = — II p(x)\ogp(x)dx 

J J\x\>l 



d6 
J\ 



^± p -^Xo g £-± p -^)pdp 

lIT Z7T 



a + l + (a-l)log^_ 



a - 1 



(12) 



p(x) 



(a-l)\x\- a - 
2n{L 1 - a - 1) 



-> 0. 



(15) 



Therefore, we know that the visit probability of each 
point x is approaching to 0. That implies there are no 
point can be visited twice by the T random walkers. In 
another word, we obtain an even distribution on those T 
visited points. So the Shannon entropy can be calculated 
by the following formula: 



T 

s = -E^ lo g^ = lo g T 



(16) 



Instead of formula EqfT^l In this case, we can estimate 
the total number of walkers (or total time steps of one 
random walker) as follows. 

We know for each "jump out and return" cycle, the 
average distance traveled by the walker is: 



D = 



l dx 



1 



L>|5|>1 

dO / p-p 

2ix r-i r 2-a 



^pdp 



(17) 



^-(1-L 1 -") 2- a 



L, 



And we know the total constraint distance can be written 
as. 

W = cL, (18) 
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where, c is a constant. In the main text, we set c = 1 or 10 
in the simulations. So, the average total number of time 
steps is: 



T= W = c(2-a) 
D l-a ' 



(19) 



Bring it into Eq |161 we get: 



S(a) = log T = log 



c(2 - a) 
l-a 



(20) 



And we know: 



dS(a) _ l-a (a- 2)c 
da ~ {2-a)c\a- l) 2 



1 — Q 



(21) 



Which is always larger than when < a < 1. That 
means S(a) is a monotonic increasing function. So, sum- 
marizing cases i and ii, we know, 



dS(a 
da 



< if a > 1 
> if a < 1 



(22) 



That implies S(a) can get its maximal value when a — 
1. Therefore, we have proved that when return probability 
r = 1 , the visit entropy S can get its maximum when a = 1 
which is consistent with our simulation result. 

Never Return (r = 0). In this case, the random walker 
will keep jumping away from the home site and never come 
back again until its energy is consumed up. 

As in the case of r = 1, assume that in each time t the 
visit times of a given site a; is a random number n(x, t) 
which follows — 1 distribution, that is: 



n(x, t) 



1 with probability p{x, t) 

with probability l—p{x,t), 



(23) 



where, p(x, t) is the visit probability of the given site x 
at the t th time step. It can be viewed as an independent 
stochastic process for given x. (We will get its expression 
latter). Then after T time steps, this site x will be visited 
£(x) times, 



(24) 



t=i 



So, its average value is: 



«(2)> = = !>(*>*)• 



(25) 



Then the Shannon entropy can be calculated as: 



?OT^ log( ^ )) 



log E 



LxL lt:(-'\\ LxL 



(26) 



= -£{ 



LxL T 



t=0 

LxL ,^ ,s T 

= - E { bg e p(s , « ) > + log r 



LxL T 



t=0 
T 



= ~f J2 {J2 p(x,t) log J2p(x,t)} + \ogT 

x t=0 t=0 

Finally, we get a concise expression: 

LxL 

T 



^ LxL 

S = - - E *) log *) + log T 



Where, 



(27) 



(28) 



t=i 



Therefore, the Shannon entropy is just the time average 
value of the entropy of the visit probability plus a 
constant logT. So, S is the function of T. However, as we 
know, T and p(x, t) are the functions of the exponent a. 
Thus, we should find the expressions T(a) and fj,(x,t,a) 
to solve the problem. However, it is very hard to get the 
mathematical explicit expressions. So we will only give 
the numeric results instead. 

i. T{a) 

For given a and W (the total distance constraint), T 
is a random variable. We know that in each time step, 
the random walker will jump out a distance l t which is 
a random number with the distance distribution pi t (I) oc 
l~ a , so according to the distance constraint, we should 
have the following inequality: 



T+l 



E h < W and E h > W - 



(29) 



t=i 



t=i 



t=i 



Then, the problem becomes a classic problem of the first 
passage time of one-sided Levy flight. We can convert 
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Fig. 6: The Relationship between W and T in different a. The 
slopes (7) of the straight lines show systematic changes 



While, the exponent 7 is not always a — 1 but changes 
with a in a relationship as Fig. [7] shows. 
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this problem into a following equivalent one: on a one- 
dimensional line, a random walker starts from the origin 
and performs the one-sided Levy motion (which means the 
random walker can only move to right but never back) 
until it will be attracted by a wall which locates at W. So 
the time T is just the first passage time to the location 
W. 

This problem has been fully discussed in the references 
[23 , 24] , and we can obtain the distribution of the random 
variable T . While, here, we only need to use the average 
value of T as the approximation. We know that [23| [24] . 



(T) oc W a -\ for l<a<2. (30) 

However, in our case a can be larger than 2. Thus, 
we cannot use this analytic result. Instead, we have done 
a large number of simulations about this one sided Levy 
flight random walk and found that the power law relation- 
ship between T and W is always hold (see also FigJBJ): 



(T) oc W~<. 



(31) 



We can use a binomial equation to approach the simula- 
tion results, so we get the following relationship: 



T(a) ( X i^ - 124106 + - 790568 ( Q ~ 1 )~ - 199708 ( a - 1 ) 2 
for 1 < a < 3, 



(32) 



Where, T is just the traveling time under the distance 
constraint W that we will use. 



Fig. 7: The Relationship between 7 and a — 1 in the one-sided 
Levy flight simulations 



after the Fourier transformation, so the solution can be 
written as 1251: 



exp[i{k -x) - |fc| Q "4]dfc. 



(34) 



It is the 2-dimensional Levy stable distribution. We know 
that it has no analytic solution therefore we will only give 
its numeric results. Finally, 



H(x, t) 



T 

ff ex P [*(* ' f ) _ \k\ a -H]dk. (35) 



4tt 2 



t=i 



Then, bring this result into the Eq|27l we can get the 
numeric result of dependency of S on a shown in Fig. |8l 

In Figil we set L = 1000, W = 10, 15, 20L. We can 
observe that the curve can get its peak at a = 1.8 which 
is close to the simulation result. Although the shape of 
the curve is different from the simulation result (Fig. 3) 
because lots of approximations are adopted in this analy- 
sis, their main features are similar. We can know that as 
the W and L increase, the optimal exponent will approach 
to the simulation result. 
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